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A  radiation  conversion  and  radiation  transport  model  applicable  to  gadolinium-exposure 
neutron  Imaging  Is  formulated.  The  resultant  linear  optical  density  representation  is  found 
to  be  In  good  agreement  with. experiment.  The  distinct  conversion-electron  contributions 
thus  Identified  confirm  the  supposition  aboutlhc  dominant  mode  of  Image  formation  in 
gadolinium-exposure  neutron  imaging;  : 


INTRODUCTION 

The  increasing  availability  ol  intense  neutron  sources 
has  contributed  to  a  broad  interest  in  experimental  neu¬ 
tron  radiography.  Recently,  studies  have  been  re¬ 
ported  wh'ch  elucidate  the  mathematical-geometrical  de¬ 
scription  of  the  optical  Image. w  These  studies  suggest 
that  a  need  exists  for  a  model  of  the  dominant  radiation 
conversion  and  transport  processes  which  contribute  to 
the  optical  film  density.  Herein  we  formulate  and  test  a 
simplified  radiation  conversion  and  radiation  transport 
model  specifically  relevant  to  gadolinium-exposure  neu¬ 
tron  radiography. 

The  experimental  arrangement  of  the  neutron  radio¬ 
graphic  system  considered  ts  schematically  represented 
in  Fig.  1.  Neutrons  are  assumed  to  pass  largely  unat¬ 
tended  through  the  film  emulsion  and  aluminum  foil 
and  thereupon  become  strongly  attenuated  in  the  gado¬ 
linium  conversion  foil.  The  neutron-nucleus  absorption 
process  in  the  converter  foil  thereupon  produces  a 
source  of  conversion  radiation  which  subsequently  leads 


to  film  blackening.  The  discussion  on  the  type  of  radia¬ 
tion  which  eventually  leads  to  the  optical  density  in  the 
film  is  deferred  to  the  end  of  this  paper. 

RADIATION  TRANSPORT  MODEL 

Ws  consider  the  energy-dependent  neutron-nucleus  in¬ 
teraction  density  at  the  point  P  in  the  converter.  If  we 
designate  this  quantity  by  F( x,),  then,  In  this  one-di¬ 
mensional  representation,  tre  may  write 

jF(x,)= j’X(E)<p(Eyx,)dE 

=  f’  SteMF.x^expt-SdZ)  Or, -*,)]<!£ 

=  £<f>(xl)exp£-.£(x1-xj)].  (1) 

In  this  expression  we  use  £  as  the  appropriately  energy- 
averaged  macroscopic  absorption  cross  section  for  neu¬ 
trons  in  gadolinium  foil  and  <*>  (x,)  represents  the  scalar 
neutron  flux  on  the  xt  plane. 

To  obtain  an  expression  for  the  production  density  of 
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FIG.  1.  Geometric  and  material  arrangement  used  for  neutron 
radiographic  imaging. 


conversion  radiation,  designated  by  R(xt),  we  assume 
that  tliis  quantity  is  directly  proportional  to  the  neutron- 
nucleus  interaction  density: 

R.(xt)=  l^Ffri^yjS^CeiJexpl  -SG^  —at,)].  (2) 

In  this  equation  we  introduce  the  subscript  "j"  to  empha¬ 
size  that  distinct  conversion  radiation  energy  groups 
may  be  produced.  Thus,  Y,  is  equal  to  the  yield  frac¬ 
tion  of  the  jth  radiation  energy  group  produced  in  this 
conversion  process. 

The  radiation  flux  at  xt  which  contributes  to  the  optical 
density  in  the  film  emulsion,  defined  by  J,(x2),  will  de¬ 
pend  upon  the  transport  kernel  for  the  specific  radiation 
energy  group  in  question.  If  we  symbolically  define  this 
positive  dimensionless  function  by  K,  ix,~'x2~‘x2),  we 
may  combine  this  kernel  with  the  radiation  production 
density  rate  by  an  integration  over  the  gadolinium  thick¬ 
ness  to  yield  an  expression  for  the  radiation  flux  at  x2: 

^(xa)=jf'i  Y}  S^(xj)exp[— S(x, -x,)]Ki(x,-xl^xa)dx1 

=  Yj  2tp(xt)f’lexp[-E(xi -xt)]ICi{x,~xt',,xa)dxi. 

*'  (3) 

In  this  expression  it  is  necessary  to  emphasize  the  lim¬ 
ited  range  of  permissible  integration  of  the  transport 
kernel  for  the  ;th  radiation  energy  group;  for  this  reason 
we  specify  the  lower  limit  of  integration  by  xr 

In  the  present  analysis  we  assume  that  the  kernel  for 
our  radiation  transport  model  can  be  adequately  approx¬ 
imated by  a  positive  definite  linear  relationship  of  the 
form 

Kj(xt~x{~x2) 

=  |l  *• ~ *|)l  !1~(pAl/«Ai.P0r,-x,)!,(4) 

whcie  RatJ  and  RAij  are  the  ranges  in  gadolinium  and 
aluminum,  respectively,  in  units  of  mass  per  area;  the 
densities  of  gadolinium  and  alvmimtn  are  represented  by 
Pod  an^  Pai  • 

If  the  radiation  groups  can  be  considered  to  be  contained 
in  a  relatively  narrow  energy  range,  then  we  may  use 
the  common  model  that  the  range,  when  expressed  in 
units  of  mass  per  area,  is  independent  of.the  material. 
Hence  we  write 


I 

I 

I 

Ra<i,j  =  b\i,)=RJ-  (5)  | 

Upon  substitution  of  Eq.  (5)  into  Eq.  (4)  and  subsequent¬ 
ly  into  Eq,  (3)  we  obtain  an  explicit  expression  for  the 
radiation  flux  at  x,: 

J,b2)=  Yj  Z<f>(xi)  f** exp[- Z(x2 -*,)]  j 

x  J  ^ri.-Pai(x2-*i)  J  -  j 

=  [i,i<^Oft)exp(-Sx1)/RjK^-pAlCv2-x1)]  ’•  j 

I 

x  -  Po^^fexptSxj)  -  expfSx,)]  \ 

f 

+  Pode*P(S*i)(*t  - 1/ Z)-p0iexp(Zx1)(xl  - l/s)}.  ! 

(6)  ; 

This  expression  for  the  radiation  liux  at  x2  may  be  sym¬ 
bolically  written  as  i 

(7)  1 

where  G(s,x1,x/,p0d)  is  a  function  of  the  gadolinium  . 

converter  only.  Since  the  radiation  flux  at  x2  will  sub-  ■ 
sequently  impinge  on  the  film  emulsion  then-assuming 
that  film  blackening  occurs  in  the  linear  exposure  \ 

range— the  resultant  optical  density  will  be  proportional 
to  the  radiation  flux.  Therefore,  the  experimentally 
measured  optical  density  resulting  fvom  the  jth  radia-  ; 
tion  energy  group,  designated  by  D.,  can  be  represented  j 
by 

Ri~ijYj\.(i>l}ci)/R*\G(Z,x2txl,paJ[R1  -pAt(x2  — Xj)],  (8) 
where  ef  is  the  emulsion  response  parameter. 

From  Eq.  (8)  it  is  evident  that,  for  a  given  neutron  ra¬ 
diographic  imaging  system,  the  optical  density  Df  for 
the  jth  radiation  energy  group  will  decrease  linearly 
with  increasing  alunimm  thickness.  If  there  are  several 
such  energy  groups,  then  the  resultant  optical  density 
constitutes  the  summation  effect  of  such  linear  contri¬ 
butions.  Thus,  the  individual  radiation  energy  groups 


FIG.  2.  Relative  optical  film  density  as  a  funciloo  of 
aluminum  thickness.' 
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TABL£  1.  Range  and  optical  density  contribution 
of  conversion  radiation. 


Radiation  Conversion  Relative  contribu-  Optical  density 
energy  radiation  tlon  to  optical  contribution 
group  range,  A  -  density,!^  (normalized  to 

O')  (mg/cm3i  (arbitrary  scale)  100%) 


5. 5  ±0,3 
18±2 
32±8 


0,599±  0.025 
0. 102  ±0.020 
0.073  ±0.016 


should  be  experimentally  dlscernable  If  they  are  suffi¬ 
ciently  distinct. 

EXPERIMENTAL  RADIATION  IMAGING 

To  test  the  validity  of  the  analysis  leading  to  Eq.  (8),  wo 
conducted  a  series  of  neutron  radiographic  experiments 
using  the  imaging  system  schematically  represented  in 
Fig.  1.  The  converter  consisted  of  a  25-/±-thick  high- 
purity  gadolinitm  foil.  Aluminum  foils  were  assembled 
to  vary  in  thickness  from  2  to  50  mg/cm*.  A  neutron 
flux  of  -5X10*  /cma/sec  was  extracted  from  the 
McMaster  University  Nuclear  Reactor  through  a  2. 5- 
em-diam  vertical  access  tube. 

Optical  film  density  is  not  an  ideal  indicator  of  radiation 
energy  group  exposure.  To  circumvent  some  of  the  po¬ 
tential  problems,  we  undertook  several  precautionary 
steps.  The  effects  of  developer  conditions  were  mini¬ 
mized  by  batch  development.  From  preceding  exposure 
calibration  measurements,  it  was  possible  to  select  ir- 
radlat’on  times  which  ensured  that  the  resulting  film 
blackening  would  fall  within  a  narrowiy  defined  range. 
The  ratios  of  times  required  to  obtain  the  same  film 
blackening  as  that  required  for  a  standard  was  taken  as 
a  relative  measure  of  optical  film  density.  The  useful¬ 
ness  of  this  latter  step  rests  on  the  assumption  that 
"reciprocity  failure”  does  not  exist  for  the  image  for¬ 
mation  process;  this  has  been  found  to  be  generally  ap¬ 
propriate  for  neutron  radiography.  * 

The  experimental  results  from  the  neutron  radiographic 
measurements  are  shown  in  Fig.  2.  It  is  graphically 
evident  that  the  relative  optical  density  Is  resolvable  In¬ 
to  three  distinct  linear  components.  By  the  methods  of 
least  squares,  we  have  calculated  the  effective  ranges. 
Rs,  of  these  three  radiation  energy  groups  together  with 
their  relative  intensities  at  xy  These  results  are  sum¬ 
marized  in  Table  I, 

Equation  (8)  can  be  reduced  to  yield  Information  on  the 
product  of  the  emulsion  performance  parameter  and  the 
conversion  radiation  yield  given  by  t}Yt.  For  this  calcu¬ 
lation  we  specify  zero  a’onlmvn  thickness  and  impose  a 
lower  limit  of  Integration  by  setting 

x^ -sxj — Rj/pQ4  ,  (9) 

Equation  (8)  now  reduces  to 

i  -id  -c-M/ttj]},  (10) 

where 

aj~R;T$/pQtl  (-1) 

with  Rj  satisfying 

RJ<(y1-x0)pQ4.  (12) 


Rearranging  Eq.  (10)  yields  an  explicit  expression  for 
the  parameter  (,Y}: 

(13) 

Herein  we  have  normalized  the  optical  film  density  to 

The  calculations  leading  to  the  results  enumerated  in 
Table  I  have  provided  both  and  the  values  used 
for  S  and  p0d  were  1400  cm"*  and  7. 8d  g/cm*.  Using 
these  parameters,  we  have  calculated  the  parameter 
tjYj,  see  Table  n. 

DISCUSSION  AND  CONCLUSION 

The  system  and  procedures  used  in  this  investigate  i 
are  typical  of  those  used  in  direct-jxposure  neutrr.n  _a- 
diography  when  gadollnivin  Is  used  as  the  converter  me¬ 
dium.  The  experimental  results,  therefore,  provide  ev¬ 
idence  concerning  the  radiation  transport  and  radiation 
recording  process  which  characterize  such  .adiographlc 
techniques. 

It  is,  however,  informative  to  compare  our  experimen¬ 
tal  result  obtained  from  this  work  with  a  result  from  a 
recent  beta-spectroscopy  study  which  employed  gadoltn- 
'  -•n  exposed  to  a  neutron  beam. 10  From  this  investiga- 
■  .on  by  Feigl  and  Rauch  it  Is  also  clear  that  three  con¬ 
version-electron  groups  dominate  the  radiation  spec- 
tnan.  We  have  extracted  the  intensities  of  the  lowest 
two  energy  groups  for  the  several  gadolinium-foil  thick¬ 
nesses  and  find  that  their  ratio  varies  from  6. 0  to  8. 0 
with  an  average  of  6. 8. u  The  most  appropriate  data 
from  «ir  work  corresponding  to  this  ratio  are  given  by 
the  performance-yield  parameters,  Table  II: 

(14) 

The  agreement  between  these  indeoendent  results  sug- 
gests  that  the  emulsion  performance  parameter,  (jt 
may  not  be  strongly  energy  dependent  for  the  range  of 
energies  involved  and  that  the  conversion  electrons  rep¬ 
resent  the  dominant  form  of  radiation  which  contributes 
to  film  blackening  in  gadolinium-exposure  neutron 
radiography. 

Based  on  the  conclusion  that  the  conversion  electrons 
represent  the  dominant  radiation  in  this  neutron  radio- 
graphic  process,  one  can  use  the  experimentally  mea¬ 
sured  radiation  ranges,  Table  I,  to  calculate  the  ener¬ 
gies  of  the  electron  groups.  We  have  evaluated  these 
energies  using  the  Katz-Penfoldu  empirical  correlation 
and  find  that  the  energies,  in  order  of  decreasing  optl- 


TABLE  U.  Performance-yield  parameter  for  each  conversion 
radiation  energy  group. 


iatlon  energy  group 


Performance-yield 
parameter,  c 
(normalized  to  100%) 

83.2% 

11.8% 
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cal  density  contribution,  aiv  give’  by  64,  135,  and 
190  keV. 

Finally,  we  point  out  that  the  design  of  improved  imaging 
devices  for  gadolinium-expos-  re  neutron  radiography 
should  be  guided  by  the  experimental  results  found  i  x 
this  investigation.  The  parameters  which  appear  to  be  of 
greatest  importance  are  the  various  contributions  to  op¬ 
tical  density  of  the  several  radiation  energy  groups, 
i.  e. ,  71.8%,  19.4%,  and  8.8%,  together  with  their 
corresponding  ranges  given  by  5. 5,  18,  and  32  mg/cm8. 
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